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We performed inelastic neutron scattering measurements on single crystals of NdFe3( 11 BOs)4 
to explore the magnetic excitations, to establish the underlying Hamiltonian, and to reveal the 
detailed nature of hybridization between the 4/ and 3 d magnetism. The observed spectra exhibiting 
a couple of key features, i.e., anti-crossing of Nd- and Fe-excitations and anisotropy gap at the 
antiferromagnetic zone center, are explained by the magnetic model including spin interaction in 
the framework of weakly-coupled Fe 3+ chains, interaction between the Fe 3+ and Nd 3+ moments, 
and single-ion anisotropy derived from Nd 3+ crystal field. The combination of the measurements 
and calculations reveals that the hybridization between 4/ and 3d magnetism propagates the local 
magnetic anisotropy of the Nd 3+ moment to the Fe 3+ network, leading to the determination of 
the bulk structure of both electric polarization and magnetic moment in the multiferroics of the 
spin-dependent metal-ligand hybridization type. 

PACS numbers: 75.10.Dg, 75.25.-j, 75.85.+t 

I. INTRODUCTION 

Coexistence of magnetic order and electric polar¬ 
ization, multiferroicity , has become a major topic 
over the past decade in condensed matter physics. 

Since multiferroicity was originally discovered in per- 
ovskite TbMnOaji various multiferroic compounds 
have been found, including J?Mn0 3 ( R — Eu, 

Gd, Tb, and Dy ),- 2 Bao.sSri.sZ^Fe^C^,” M 3 V 2 O 8 P 
CoC^Cbj— MnWC> 4 )£ CuFeCb^ LiCu 202 <~ LiCuVCbf— 
and Ba 2 CoGe 2 C> 7 *i£ Recent theoretical and experimen¬ 
tal studies revealed that the electric polarization in these 
compounds is driven by magnetic long-range order .Air— 

Since the structure of the order is determined by the 
exchange pathways and the magnetic anisotropy, experi¬ 
mental identification of the magnetic Hamiltonian is very 
important for understanding multiferroics. 

The rare-earth ferroborates i?Fe 3 (B 03)4 (R = rare- 
earth metal) are a series of new multiferroic compounds 
containing R 3+ (4 f n ) and Fe 3+ (3d 5 S = 5/2) as mag¬ 
netic ions. The variety of the magnetic anisotropy of the 
R 3+ moments ( R = Y, Pr, Nd, Sm, Gd and Tb) com¬ 
bined with the interaction between the Fe 3+ and R 3+ 
moments (f-d coupling) gives rise to diverse magnetoelec¬ 
tric (ME) effects as a function of the R 3+ ions— r— In 
these compounds the mechanism of magnetoelectricity is 
explained by the spin-dependent metal-ligand hybridiza¬ 
tion model : 22 : 23 

The crystal structure has the trigonal space group 
i?32, which belongs to the structural type of the min¬ 
eral huntite CaMg 3 (C 03)4 — As shown in Fig. Oja) the 
main feature is that distorted FeOg octahedra form spi¬ 
ral chains with threefold screw-axis symmetry along the 


crystallographic c - axis. Each chain includes three Fe 3+ 
ions in the unit cell. The chains are separated by the 
R 3+ and B 3+ ions. 

In NdFe 3 (BC> 3)4 the Nd 3+ ions (4/ 3 ) carry magnetic 
moment with J = 9/2. The magnetic susceptibility 
showed anisotropic decrease below 29 K, and the heat 
capacity showed well-defined A type anomaly at the same 
temperature, implying a phase transition to an antifer¬ 
romagnetic (AF) ordered state with Neel temperature of 
T/v = 29K— AtT> Tn the susceptibility is concave 
downward, indicating the short-range AF order because 
of the low dimensionality of the magnetic system. Spon¬ 
taneous electric polarization simultaneously appears in 
the AF ordered phased The electric polarization signif¬ 
icantly increases upon applying a magnetic field parallel 
to the a - axis. The magnitude of the electric polarization 
reaches P a ~ 300 pC/m 2 at 1.3 T and 4.2 Kj 17 ' 18 which 
means that the magnetization along the a - axis induces 
the large electric polarization along the a - axis. A neu¬ 
tron diffraction study exhibited an easy-plane type AF 
order at T <Tjq-, the Fe 3+ and Nd 3+ magnetic moments 
align ferromagnetically along the a - axis and propagate 
antiferromagnetically along the c - axis with the propa¬ 
gation vector k = (0, 0, 3/2) in Fig. [Qa) i 26 i 27 Both of 
Fe 3+ and Nd 3+ moments are simultaneously ordered at 
the Tn, indicating non-negligible f-d coupling. Further 
decreasing the temperature at T < Tic = 13-5 K, the 
commensurate (C) magnetic peak splits into a pair of in¬ 
commensurate (IC) peaks where the magnetic moments 
are in the ab - plane and the AF helix propagates along 
the c - axis. 

Magnetic dynamics have been investigated by spectro¬ 
scopic methods using electromagnetic waves. ESR mea- 
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FIG. 1: (Color online) (a) The crystal structure and magnetic 
structure of NdFe 3 (BC> 3)4 (hexagonal, space-group R32). (b) 
Laue image of a crystal with the c* - axis parallel to the 
incident X-ray beam. Crystals of NdFe3( 11 B0a)4 on the alu¬ 
mina holder viewed from above (c) and from the side (d). 
(e) Reciprocal space in the a* - c* plane. The square at 
( h,k,l ) = (— 1, 0, — 1) is a nuclear Bragg reflection. The cir¬ 
cle at (0,0,—1.5) is a magnetic Bragg reflection. Red curves 
indicate spectra measured with fc;/ /a* for hui = 0,1 and 3 
meV. The blue shaded area indicates the range of the ob¬ 
served scattering plane for hi = 2.352 A -1 and tku = 0 meV 
when sample is rotated. 


surement detected an energy gap suggesting a uniaxial 
magnetic anisotropy in the ab - plane .- 28 It also detected 
lifting of the Kramers doublet of the Nd 3+ ion due to the 
molecular field from the neighboring Fe 3 " 1 " ions. Opti¬ 
cal spectroscopy provided the energy levels of the crystal 
field of the Nd 3+ ion and determined the parameters of 
the crystal field Hamiltonian^ 

In NdFe 3 ( 1103)4 exhibiting the strong f-d coupling, 
the investigation of excitation spectra including Fe 3+ spin 



FIG. 2: (Color online) Inelastic neutron spectra projected 
onto c* - axis at (a) 41 K and (b) 15 K. The incident neutron 
energy was 11.46 meV. 


wave and Nd 3+ crystal field in a wide wave-vector en¬ 
ergy space is crucial in order to identify the magnetic 
Hamiltonian and to unravel the detailed nature of hy¬ 
bridization between 3 d and 4/ magnetism. Furthermore 
in constructing the Hamiltonian, careful consideration 
of magnetic anisotropy is important in the nrultiferroics 
of the spin-dependent metal-ligand hybridization mecha¬ 
nism typef 14,1 ° in which the magnetic anisotropy directly 
determines the polarization structure. 

In the present paper we study inelastic neutron scatter¬ 
ing (INS) spectra on NdFe 3 ( 11 BOs )4 to explore the mag¬ 
netic excitations and to establish the underlying Hamil¬ 
tonian. Following to the introduction we describe the 
experimental details about the sample preparation and 
the setup of INS measurements in Sec. |H] Subsequently 
in Sec.|TTT]the INS spectra of NdFe 3 ( 11 BOs )4 are demon¬ 
strated. We observed spin waves of the Fe 3+ moment 
below 6 meV and transition between the lifted states of 
Kramers doublet of the Nd 3+ ion at 1 meV. A couple 
of characteristic features are an anti-crossing of the Fe- 
and Nd-excitations, and a small anisotropy gap at the 
antiferromagnetic zone center. In Sec. IIVI the magnetic 
model including an in-plane anisotropy derived from the 
crystal field excitation of the Nd 3+ moment and the non- 
negligible f-d coupling is constructed. The observed 
spectra are successfully analyzed by the linear spin wave 
theory based on the model. The origin of the in-plane 
anisotropy is revealed to be the crystal field of the Nd 3+ 
ion. In Sec. [V] possibility of magnetic anisotropy of the 
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FIG. 3: (Color online) Constant energy cuts at 0 meV (a), 1 meV (b), 2 meV (c) and 3 meV (d) in the a* - c* plane. The 
horizontal and vertical axes represent the a* and c* axes, respectively 


Fe 3+ moment is discussed. It is turned out that the 
anisotropy is very small, and, instead, the Fe 3+ mo¬ 
ment inherits an in-plane anisotropy through hybridiza¬ 
tion with the Nd 3+ moment. The conclusions are given 
in Sec. eh The magnetic Hamiltonian in NdFe 3 ( n B 03)4 
is established in the present study. Combination of the 
measurement and the detailed calculation revealed that 
the hybridization between 4/ and 3 d magnetism prop¬ 
agates the local magnetic anisotropy of the Nd 3+ ion 
to the Fe 3+ network, resulting in the bulk structure 
of multiferroics. The local symmetry of the rare-earth 
ion is a driving force for the non-local multiferroicity in 
NdFe 3 ( 11 B0 3 ) 4 . 


II. EXPERIMENTAL DETAILS 

Single crystals of NdFe 3 ( n B 03) 4 were grown by a flux 
method^ We first synthesized polycrystalline samples 
from the starting materials, Nd 203 , Fe 2 C> 3 , and 11 B 2 03 . 
The stoichiometric amounts of the starting materials 
with a total mass of about 16 g were mixed, ground, and 
put into an alumina crucible. The crucible was heated 
at 980 °C for 72 h. The flux is Bi 2 Mo 3 Oi 2 + 3 n B 2 0 3 
+ 3/5 INchOs; Bi 2 Mo 30 i 2 was synthesized by the solid 
state reaction from B^Os and M 0 O 3 inside an alumina 
crucible at 600 °C for 24 h. A mixture of about 60 g 


of NdFe 3 ( n BC> 3) 4 and the flux with the mass ratio of 1 
: 3 was put into a platinum crucible inside the alumina 
crucible. The crucible was heated to 1000 °C for 4 h, 
kept at this temperature for 1 h, cooled to 962 °C for 1 
h, and slowly cooled down to 870 °C for 120 h; then the 
furnace was shut down to the room temperature. The 
flux was removed by decanting at 900 °C, and washing 
the crystals with HC1 solutions. 

We coaligned 22 pieces of single crystals so that the 
crystallographic a* - c* plane is horizontal. Alignment 
was performed by transmission Laue method using a 
high energy X-ray Laue camera. The X-ray source was 
YXLON MG452 and the maximum energy of the white 
X-ray beam was 310 keV. We recorded Laue patterns us¬ 
ing a high-speed CCD camera, with imaging size 10 cm 
x 10 cm (1024 x 1024 pixel). Figure [Ub) shows a Laue 
image of a crystal in the c* - plane. This pattern exhibits 
the threefold symmetry along the c* - axis. We placed 
the crystals on an alumina holder as shown in Figs.[ljc) 
and Ed). The average mass of the crystals was 0.1 g. 
The total mass of the sample was 2.1 g. 

The INS experiment was performed at the High Res¬ 
olution Chopper Spectrometer (HRC) installed in the 
Material and Life Science Experimental Facility of J- 
PARCi^ic — 3 At the HRC white neutrons are monochrom- 
atized by a Fermi chopper synchronized with the produc¬ 
tion timing of the pulsed neutrons. The energy transfer 
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hid was determined from the time of flight (TOF) of the 
scattered neutrons detected at position sensitive detec¬ 
tors (PSDs). The To chopper was set at 50 Hz, a colli¬ 
mator of 1.5° was installed in front of the sample, and the 
“S” Fermi chopper with 200 Hz was used to obtain high 
neutron flux. We used a GM-type closed cycle cryostat 
to achieve 41 K and 15 K. The energy of the incident 
neutron beam was £) = 11.46 meV yielding an energy 
resolution of A E = 0.3 meV at the elastic position. 

Figure QJe) illustrates the a* - c* scattering plane. 
Reciprocal lattice positions at q = (—1,0,—1) and 
(0,0,—1.5) are the positions of nuclear and magnetic 
Bragg peaks, respectively. Throughout this paper q is 
expressed in reciprocal lattice unit, q = ( h,k,l ). INS 
spectra with kij/a* were measured at 41 K and 15 K. 
Red curves in Fig. []Je) indicate the measured q-ranges 
for Hu = 0,1 and 3 meV. At 15 K INS spectra that cover 
wide q-range were measured by rotating the crystal by 
70 degree in 2 degree steps. The q-range in the scan 
for hid = 0 meV is indicated by the blue shaded area 
in Fig. [He). The range of out of plane momentum is 
|q_i_| < 0.41 A -1 . In the following, the out of plane spec¬ 
tra are integrated in the central range |qj_| < 0.13 A -1 , 
and we show all spectra in the (ft., 0, Z) plane. 


III. EXPERIMENTAL RESULTS 

The INS spectra projected onto the c* - axis at 41 K 
and 15 K are shown in Figs.[U(a) and[2](b), respectively. 
At 41 K a diffuse spectrum of paramagnetic scattering 
emerges from l = —1.5. At 15 K dispersive excitations 
emerge in the energy range of 2.5 meV < hid < 4.5 meV, 
which is interpreted as spin waves of the Fe 3+ moments. 

Figures[3Ka)-[3Kd) display neutron spectra at 15 K sliced 
at the energies of 0, 1, 2, and 3 meV in the a* - c* plane. 
The white arcs are because of the absence of the neutron 
detectors between the detector banks. In Fig. (Ha) the 
peak at q = (—1,0,-1) is a nuclear Bragg reflection, and 
the peak at (0,0,—1.5) is a magnetic Bragg reflection. 
The peaks at other qs are not identified; they may be 
Bragg reflection from minor grains of crystals. The rings 
expanding from (0,0, —1.5) in Figs. E{b)-[3Jd) imply that 
a dispersive excitation appears from (0,0,—1.5). The 
rings are flattened along the c* - axis, which means that 
the dispersion along the c* - axis is steeper than that 
along the a* - axis. This is consistent with the naive 
prediction from the crystal structure that the intrachain 
interaction along the c - axis is strong compared with the 
interchain interaction. 

The INS spectrum at 15 K projected onto huj - (0,0, l) 
plane by the integrating the neutron intensity in the 
ranges of —0.1a* < q < 0 .1a* along the a* direction in 
the scattering plane and —0.17a* < q < 0.17a* perpen¬ 
dicular to the scattering plane is shown in Fig. Eta). The 
spectrum projected onto hid - (h, 0, —1.5) plane by the in¬ 
tegration in the ranges of —0.17a* < q < 0.17a* perpen¬ 
dicular to the scattering plane and —1.6c* < q < —1.4c* 
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FIG. 4: (Color online) Color maps of the inelastic neutron 
scattering spectra obtained at HRC, along the (a) l and (b) 
h directions at 15 K. The open circles are the peak positions 
extracted Gaussian fits and the curves are the calculated spin 
wave dispersions. Fitting curves reasonably reproduce the 
experimental data. White area in the middle in the panel (a) 
is the gap between neutron detectors banks. 


along the c* direction is shown in Fig. EKb). We clearly 
observed the spin waves of the Fe 3+ moments around 
q = (0,0,—1.5), and the flat excitation at about 1.0 
meV which is the transition between the lifted states of 
Kramers doublet of the Nd 3+ ion. The spin waves of the 
Fe 3+ moments are more dispersive along the c* direction 
than along the a* direction, which is consistent with the 
flattened ring in Figs. E£b)-[3Id). 

A series of hid dependence of the neutron intensity ob¬ 
tained by the integration in the same q-ranges that were 
used for the display of Figs. EK a ) and Elb) were fitted 
by Gaussian functions to investigate the detailed struc¬ 
tures of the excitations. The peak energies were plotted 
as open circles in Figs. EKa) andEKb). These data will 
be used in the analysis section. The white circles around 
q = (0, 0, —1.5) exhibit an anti-crossing between the spin 
wave of the Fe 3+ moments and the flat mode of the Nd 3+ 
moments, meaning that the Fe 3+ moments interact with 
the Nd 3+ moments. 

Figure ED a) shows the hid dependence of the neutron 
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FIG. 5: (Color online) (a) hui dependence of the neutron 
intensity that is obtained by the integration in the q-range 
of h — [—0.5,0] and l = [-1,-0.7]. The solid curve is fit 
to the data by Gaussian functions, (b) Constant-Q scan 
at (0,0,—1.5). The spectra were integrated in the range of 
h = [—0.1, 0.1] and l = [—1.6, —1.4]. The error bars are inside 
the symbols. The solid curves are fits to the data by Gaussian 
functions. The gray shaded areas are the energy resolutions: 
Gaussian functions with full width at half maximum of 0.3 
meV. 



FIG. 6 : (Color online) Exchange paths in NdFes(B 03 ) 4 . Ji 
and J 2 are the nearest and 2 nd nearest neighbor exchange 
interactions between the Fe 3+ moments. J 3 is the nearest 
neighbor exchange interaction between the Fe 3+ and Nd 3+ 
moments. The x - axis is parallel to the a - axis and the 2 - 
axis is parallel to the c - axis. The y - axis is vertical to the 
x and 2 - axis. 


intensity obtained by the integration in the g-range of 
h = [—0.5,0] and l = [-1,-0.7], where the Nd 3+ -based 
mode is not affected by the Fe 3+ -based spin-wave mode. 
The dashed curve is incoherent elastic scattering and the 
blue curve is a Gaussian fit. The red curve is the sum 
of the dashed and blue curves. The gray shaded area is 
the energy resolution. From the peak position we iden¬ 
tify the magnitude of the energy split of the Kramers 
doublet to be 0.98 meV. Figure [5][b) shows a constant- 
q scan at (0,0,—1.5), the AF zone center, that is ob¬ 
tained by integration in the g-range of h = [— 0 . 1 , 0 . 1 ] 
and l = [—1.6,—1.4]. The dashed curve is incoherent 
elastic scattering, the blue curves are Gaussian fits, and 
the red curve is the sum of the dashed and blue curves. 
The gray shaded areas indicate the energy resolutions. 
The magnitude of the energy gap at the AF zone center 
is estimated as 0.57 meV. The gap implies the existence 
of an anisotropy in the ab - plane. The excitation at 1.3 
meV is the transition between the split of Kramers dou¬ 
blet of the Nd 3+ ion. The peak energies of the Nd 3+ ion 
in Fig. [5](a) and Fig. [5](b) are obviously different. The 
energy difference is due to the hybridization between the 
Fe 3+ and Nd 3+ modes. The energy widths of the ex¬ 
citations at foul = 0.57 meV and 1.3 meV are broader 


than the experimental resolution. The energy split of 
0.98 meV and the energy gap of 0.57 meV are consis¬ 
tent with the magnetic excitations reported from ESR 
measurement^ and optical spectroscopy^ 

IV. ANALYSIS 

In order to identify the magnetic model realized in 
NdFe 3 ( n B 03 ) 4 , we consider the following Hamiltonian: 

n = - ^ JiS, ■ s, - ^ j 2 Si ■ Sj 

n.n. n.n.n. 

-Y'JaSi-Jk + Y.U cf( J A>, (1) 

n.n. k 

where the x - axis is parallel to the crystallographic a - 
axis and the 2 - axis is parallel to the c - axis. J\ and 
J 2 are the exchange interactions in the nearest and 2 nd 
neighbor paths of the Fe 3+ ions as shown in Figs. [Hla) 
and[5](b). These terms mainly determine the dispersions 
along the c* - axis and a* - axis, respectively. J 3 is the 
nearest neighbor exchange interaction between the Fe 3+ 
and Nd 3+ moments, which induces the anti-crossing be¬ 
tween the Fe 3+ and Nd 3+ modes. In Eq. |T]), positive 
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(negative) signs of the exchange parameters correspond 
to ferromagnetic (antiferromagnetic) exchange interac¬ 
tions. T-Lcf is the crystal field Hamiltonian of the Nd 3+ 
ion. 

There are five Kramers doublets in the crystal field of 
the Nd 3+ ion above Tn = 30 K. The first excited en¬ 
ergy is about 8 ineV. and second one is about 17 meV 
as reported in Ref. 29. In the calculation of the low- 
energy excitations below 6 meV it is assumed that only 
the ground state in the crystal field of the Nd 3+ ion hy¬ 
bridizes with the Fe 3+ moments. For this assumption, 
we introduce the ground state of the crystal field of the 
Nd 3+ ion. 

At the Nd 3+ site with D 3 symmetry, the crystal field 
Hamiltonian can be defined as follow: 

Hcf = BqCq + B*C 4 + *H 3 4 (C 4 + d 3 ) 

+£ 0 6 C 0 6 + iBl(Cl + C e _ 3 ) + Bl(Cl + C%). (2) 

The B y are the crystal field parameters and the C% are 
the spherical tensor operators. We used the values of the 
parameter B? reported in Ref. Matrix elements of 
J x , J v , and J z in the ground state doublet \g±) are as 
follow: 


where jj = —3.80 x 10 -5 is a Stevens’ factor. Eq. © 
means that the classical energy of the crystal field gives a 
six-fold anisotropy. Since the sign of the Bq'jj is negative, 
the easy-axis is the a - axis that is consistent with the 
magnetic structured The magnitude of the calculated 
anisotropy energy is 65.2 geV. We effectively include the 
six-fold anisotropic energy as —D (J x ) 2 in the Hamilto¬ 
nian Eq. ©■ The coefficient D is 23.5 /reV, which is 
defined by the relation: 



V231 

~t 6 ~ 


BqIJ ( £ 


x 2 = 65.2 geV. 


(7) 


Next the operators of total angular momentum J are 
approximated as the operators of pseudo-spin s = 1/2 
because the ground state is Kramers doublet and the de¬ 
gree of freedom is two. The connection between operators 
of the total angular momentum J and the pseudo-spin 
s is determined by the matrix elements Eq. (©. In this 
approximation the operators of the total angular momen¬ 
tum is redefined as 


J x =f3s x , J v =f3s y , J z =as z . (8) 


' (9t\J x \9±) = - l- 6 ? 

< (StI JV ls±) = (3) 

(g ±I J z I g±) = t0.98 = =f^, 

<s±l J x l<?±) = (9 ±I j y I g±) = (5tI J z I g±) = 0 , 

meaning that the total angular momentum is anisotropic, 
favoring in-plane. Next we inspect anisotropy within the 
ab - plane. The quantization axis is transformed from 
the 2 - axis to the x - axis because the direction of the 
spin in the order state is along the x(a ) - axis. Then, the 
matrix elements of J x , J v , and J z in the redefined state 
\g±) are calculated as follow: 

' (g’ ± \J x \g' ± )=T^ 

< <4l J!/ l4) = §, (4 ) 

<4|J*|4)=±fb 

<41 JX 14) = <41 J y 14) = (41 J z 14) = 0 . 

The operator Rf(Cg + C 3 6 ) leads to an anisotropy 
within the ab - plane. In order to quantify this, we ex¬ 
press the total angular momentum of the Nd 3+ ion in the 
ground state doublet as 

J x = — cos 9, J v = — sind, (5) 

2 2 ' 

where 9 is the angle between the a - axis and the moment. 
Thus classical energy of Hf (Cf + C- 6 ) becomes 


The Hamiltonian Eq. © is, thus, represented by 


« = -E-'iSVSj- V 

n.n. n.n.n. 

( P4 \ 

\p*l ( 0 ) 

n.n. y OLS ^ J k 

We subsequently calculate the spin wave spectrum of this 
Hamiltonian using Holstein-Primakoff (HP) transforma¬ 
tions. The HP transformations of the spin operators of 
the Fe 3+ moments and the pseudo-spin operators of the 
Nd 3+ moments are written as 



The a\. a t . b\ and bi are bosons operator in each sublat¬ 
tices. The quantization axis is parallel to the x(a) - axis. 
We introduce spatial Fourier transformation via 


Bt(C‘ + C%) = 


cos6«, (6) 


4(4 = -^^2 c i e iqVi (* = 1~24), (16) 
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where N is the number of unit cells in the system, and 
{c;} = {a i; 6 j} are the bosons operators on each sublat¬ 
tice. By using this notation we obtain 


q a 

\ Bj -j ( q H (•?)+ hx - 

ij 


(17) 


The eigenvalues of the matrix (A + B ) (A — B ) give the 
squares of the energy of the normal modes^i: 


(A + B)(A- B) Xr(q) = {M <?)} 2 Xr{q)- (18) 


V. DISCUSSION 


In i?Fe 3 (B0 3 ) 4 , the magnitude of the electric polar¬ 
ization by the R 3+ and Fe 3+ ions are locally determined 
by the magnetic moments . 22,23 In case of NdFe 3 (B0 3 ) 4 
the existence of the in-plane anisotropy favoring order 
along the a - axis by the Nd 3+ and/or Fe 3+ ions is a key 
to the emergence of the multiferroicity. In the analysis 
section, uniaxial anisotropy of the Fe 3+ moments is not 
considered for the simplicity. In this section, we discuss 
possible magnetic anisotropies in the ab - plane of the 
Fe 3+ moments. 

The conventional origin of anisotropy in Fe 3+ -based 
magnets is magnetic-dipole interaction or single-ion 
anisotropy. The magnetic-dipole interaction between 
spins Si and S 3 is represented by 


The dispersions obtained from this calculation are 
indicated by the yellow solid curves in Figs. HJa) and 
HKb). In fact we obtained 24 modes of spin waves, 
but due to the trigonal symmetry of the lattice only 8 
modes have non-zero spectral weight. The modes at 
the highest and second-highest energies around the AF 
zone center, q = (0,0,—1.5), in Fig. |U(a) are two-fold 
degenerated, and the highest energy mode in Fig. [4|b) 
is four-fold degenerated. At the zone center, there 
are two modes at Hlu = 1.12 meV and 1.23 meV, and 
another two modes at huj = 0.32 meV and 0.61 meV. 
The 0.32 meV gap is caused by the anisotropy of the 
Nd 3+ moment in the ab plane. It vanishes if the D is set 
to zero. The modes at 1.12 meV and 1.23 meV is the 
Nd 3+ level after hybridizing with the dispersive Fe 3+ 
spin waves. The small splittings of both modes are due 
to the easy-plane type anisotropy of the Nd 3+ ion, i.e., 
the effect of q / /3. These splittings are the origin of 
the observed broadenings of the experimental peaks in 
Fig. H{b). White circles in Figs.[4|a) and[4|b) are fit by 
the mean energy of the split modes. y 2 s were calculated 
for the parameter set with the step sizes of 8J\ = 0.001 
meV, <5J 2 = 0.001 meV, and S J 3 = 0.1 /xeV in the ranges 
of |Ji| < 1 meV, | J 2 I < 1 meV, and | J 3 | < 0.1 meV. The 
obtained parameters set for the minimum y 2 is listed in 
Table [I] The fit to the data provides excellent agreement 
with the overall spectrum. It should be noted that the 
anisotropy gap of about 0.57 meV at the zone center is 
quantitatively reproduced by using the fixed parameter 
of D = 23.5 /zeV obtained from reported value of the 
parameter B 51 S It is revealed that the origin of the 
in-plane anisotropy is the crystal field of the Nd 3+ ion. 


H dip = Y, & ■ S i ~ 3 ( S * • e«) & • ey)} , 

i.j r b' 

(19) 

where r,; ? and e l3 are respectively the distance and the 
unit vector along the bond between i and j. For a 
collinear in the ab - plane, the dipole-interaction energy is 
independent of the angle to the a - axis. This is due to the 
threefold screw-axis symmetry along the c - axis, which 
also dictates that further neighbor interactions vanish. 
Thus, the magnetic-dipole interaction is not the origin of 
the anisotropy in the ab - plane. 

Next, we consider the single-ion anisotropy of the Fe 3+ 
moments. There is only one inequivalent site for the 
Fe 3+ ion and the local anisotropy is uniquely determined. 
Since the screw axis 3i or 32 is along the Fe 3+ chain, the 
FeOg octahedra are transformed by 2it/3 rotation around 
the c - axis one another. Therefore, three local coordi¬ 
nates, { Xi , Y u Zi }, can be defined on the FeC>6 octahedra 
for the anisotropy. Here {z = 1, 2, 3} are the labels of the 
Fe 3+ sites. Since the Fe 3+ moments are collinear in the ab 
- plane, we discuss the anisotropy only in the ab - plane. 
Then, the general single-ion anisotropy to 4th order in 
the Fe 3+ spin-operators is expressed by 


n a 


E 


+a x i (S 

+a x 3 y ^ 


a x 2 (Sf *) 2 

+ a y 4 ( S i') 


+ a y 2 + a xy Sf' S j V * 

4 + a x 2 y 2 (E‘) 2 (E) 2 



( 20 ) 


TABLE I: Parameters obtained by the linear spin wave cal¬ 
culations 


Ji [meV] J 2 [meV] J 3 [/reV] 

D [/zeV] (fixed) 

,x 2 

-0.482 -0.054 7.9 

23.5 

0.993 


The a x 2, a y 2, a xy , a x i, a y 4, a x 2 y 2, a x 3 y , and a xy 3 are 
independent coefficients. Here we define the the local 
coordinate {Aj, Y t , Z \} as the same as the global one 
{x,y,z}. Then the relation between the spin operators 
defined on the local coordinates {A 2j3 , 1 / 2 , 3 , ^ 2 , 3 } and the 











those defined on the global one is as follows: 

| s 2 *- = -isf + ^s», 

I qU _ qx 1 qy 

* 2 - 2 2 2 2 ’ 

l Sf = Si, 

{ qX 3 _ ^ qx VZ qy 

^ “ 2* 3 T^ 3 ’ 

cW _ qx _ 1 qy 

J 3 ~ 2 ° 3 2° 3 ’ 

S 3 3 = S z 3 , 

where the x - axis is parallel to the crystallographic a - 
axis and the 0 - axis is parallel to the c - axis. These rela¬ 
tions are substituted into Eq. (l20l) to express the single¬ 
ion anisotropy in global coordinates. Hereafter we classi¬ 
cally calculate the anisotropy energy of the collinear AF 
structure in the ab - plane. The Fe 3+ spin operators are 
classically expressed by 

{Sf,Sf,SZ) = (S cos 0,S sin 0,0), 

(Si, Si, Si) = (-S cos 6 ,—S sin 9, 0), (23) 

(Si, Si, Si) = (Scos9, S B m0,0), 


( 21 ) 


( 22 ) 


disorder which breaks the threefold rotation symmetry, 
for instance, lattice distortion, and quantum and ther¬ 
mal fluctuations . 36,37 The direction of the Fe 3+ moment, 
therefore, is determined by the anisotropy of the Nd 3+ 
moment through the f-d coupling. It can be said that the 
magnetic anisotropy of the Nd 3+ moments by the crystal 
field drives the multiferroicity in NdFe 3 ( 603 ) 4 . 

The multiferroic mechanism of I?Fe 3 ( 603)4 is spin- 
dependent metal-ligand hybridization mode l 22 : 23 where 
the relation between the electric polarization and the 
magnetic moment is locally determined by the symmetry 
of O 2- ions around the magnetic ion. In a collinear mag¬ 
netic structure the local magnetic anisotropy is a cast¬ 
ing vote in the determination of the magnetic structure, 
and, consequently, in the determination of the electric 
polarization as well. In NdFe 3 ( 603)4 the crystal field 
of the Nd 3+ ion is revealed to be the origin of the mag¬ 
netic anisotropy, which determines the bulk structure of 
multiferroics. This is in contrast with the multiferroic 
materials of which the mechanism is the spin current 
modeljiii 2 where the relation is determined by the ge¬ 
ometry of neighboring magnetic moments. 


VI. CONCLUSION 


where the Fe 3+ moments are functions of the angle 9 
between the a - axis and the Fe 3+ moments in the ab - 
plane. It is found that energy is independent on the angle 
6. Consequently, the single-ion anisotropy does not give 
the anisotropy of the Fe 3+ moment in the ab - plane. 

In the multiferroic compound 6 a 2 CoGe 207 with the 
metal-ligand hybridization mechanism, it was reported 
that an interaction between the electric polarization de¬ 
termines the magnetic anisotropy in the easy-planej 35 
We, hence, consider the electric-polarization interaction 
in NdFe 3 ( 603)4 using the mechanism, where the lo¬ 
cal electric-polarization at the Fe 3+ site is expressed by 
Pi = tJ2i ( e i ■ Si) 2 ei. Here e/ is a unit vector along the 
bond to ligands (in this case oxygen), and t is a coupling 
constant related with the metal-ligand hybridization and 
the spin-orbit interaction. The polarization interaction 
is represented by JpY^ijPi ' Pj ■ Jp is ferroelectric cou¬ 
pling constant between the electric polarizations p { . The 
energy of the polarization interaction was calculated as 
a function of angle between the a - axis and the Fe 3+ 
moments in the ab - plane. It was, then, found that 
the energy is independent on the direction of the Fe 3+ 
moments in the ab - plane. Thus, the polarization inter¬ 
action does not cause the uniaxial anisotropy in the ab - 
plane. 

Accordingly, at the Fe 3+ site the conventional sources 
of magnetic anisotropy such as the magnetic-dipole in¬ 
teraction and single-ion anisotropy, and the polarization 
interaction do not lead to the a - axis anisotropy under 
the restriction that the crystal symmetry is preserved. 
This means that the Fe 3+ moment does not have any 
uniaxial anisotropy in the ab - plane unless there is any 


We performed INS measurements to explore the mag¬ 
netic excitations, to establish the underlying Hamilto¬ 
nian, and to reveal the detailed nature of hybridization 
between the 4/ and 3 d magnetism in NdFe 3 ( n 603 ) 4 . 
Overall spectra are reasonably reproduced by spin-wave 
calculation including spin interaction in the framework 
of weakly-couped Fe 3+ chains, f-d coupling, and single¬ 
ion anisotropy derived from the Nd 3+ crystal field. Hy¬ 
bridization between the 4/ and 3d magnetism is probed 
as anti-crossing of the Nd- and Fe-centered excitations. 
The anisotropy gap observed at the AF zone center is 
explained by the crystal field of the Nd 3+ ion in the 
quantitative level. Magnetic anisotropy of the Fe 3+ ion 
allowed in the present crystal structure is small so that 
it cannot be dominant. Combination of the measure¬ 
ments and calculations revealed that the hybridization 
between 4/ and 3d magnetism propagates the local mag¬ 
netic anisotropy of the Nd 3+ ion to the Fe 3+ network, 
resulting in the bulk magnetic structure. In the multi¬ 
ferroics of the spin-dependent metal-ligand hybridization 
type, the local magnetic anisotropy controls the electric 
polarization, meaning that the local symmetry of the 
rare-earth ion is a driving force for the non-local mul¬ 
tiferroicity in NdFe 3 ( 6 C> 3 ) 4 . 
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